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AN  OPTICAL  PARTICLE-SIZING  COUNTER  FOR 


IN-SITU  MEASUREMENTS 


By 

Don  Holve  and  Sidney  A.  Self 


High  Temperature  Gasdynamics  Laboratory 
Mechanical  Engineering  Department 
Stanford  University 
Stanford,  California 


ABSTRACT 

A particle  sizing  counter  is  described,  suitable  for  In-sltu  measurements 
in  two  phase  flows  of  laboratory  scale.  It  employs  near-forward'  scatter  from 
the  focus  of  a He-Ne  laser  beam,  together  with  pulse  height-analysis  of  the 
signals  from  individual  particles.  A novel  and  essential  feature  of  the 
technique  is  a numerical  inversion  scheme  to  unfold  the  dependence  of  the 
scattered  signals  on  particle  trajectory  through  the  measurement  volume.  The 
Inversion  procedure  is  performed  by  an  on-line  computer,  and  utilizes  a prior 
calibration  with  tilonodisperse  aerosols  of  known  size.— -As  presently  configured, 
the  Instrument  has  a demonstrated  capability  of  determining  size  distributions 
in  the  diameter  range  2-25  ym,  at  concentrations  up  to  cm"3.  The 

measured  dependence  of  response  on  particle  diameter  agrees  well  with  cal- 
culations from  Mle  scattering  theory.  It  is  anticipated  that  the  technique 
can  be  extended  to  cover  particle  diameters  un  to  at  least  50  ym,  and  down  to 
0.5  ym  and  concentrations  up  to  10^  cm“^.  It  should  also  be  adaptable  to  hot 
flows  and  absorbing,  irregular  particles. 
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1.0  INTRODUCTION 


The  measurement  of  particle  size  distributions  In  two  phase  flows  Is 
of  considerable  current  Interest,  especially  In  connection  with  energy  conver- 
sion devices  such  as  liquid  spray  and  pulverized  fuel  combustors,  and  parti- 
culate cleanup  devices  such  as  electrostatic  precipitators.  While  the  parti- 
culate characteristics,  such  as  mean  diameter,  size  distribution,  mass  loading 

etc.,  may  vary  quite  widely  In  such  systems,  the  sizes  of  Interest  generally 

6 “3 

fall  In  the  range  0.5  - 50  ym,  with  concentrations  up  to  10  cm  (for  the 
smaller  sizes). 

The  available  measurement  techniques  [1,2]  also  vary  widely  In  type  and 
capability.  Several,  Including  microscopy,  cascade  Impactors,  Coulter  counters, 
mobility  analyzers  and  commercial  optical  counters  require  a sample  to  be 
extracted  from  the  flow.  This  poses  problems  related  to  obtaining  a represen- 
tative sain)le,  especially  for  the  larger  sizes  and  In  hot,  high  velocity 
flows.  Also,  many  sampling  methods  are  cumbersome  and  slow  In  operation. 

For  these  reasons  optical  techniques  are  of  especial  Interest  since.  In  princi- 
ple, they  are  capable  of  making  In  situ  measurements  with  continuous  and  rapid 
readout.  Moreover,  by  using  lasers,  they  are  adaptable  to  high  temperature 
systems  having  high  thermal  radiation  background. 

Optical  techniques  all  depend  on  Mle  scattering  [3,4],  and  can  be  broadly 
divided  Into  Imaging  and  non- Imaging  types.  The  former.  Including  flash  photo- 
graphy [5]  and  holography  [6],  are  limited  In  practice  to  sizes  ~ 30  ym,  and 
pose  a difficult  data  reduction  problem.  Non-lmaglng  methods  can  be  subdivided 
Into  two  classes:  those  which  measure  on  a large  number  of  particles  simul- 
taneously, and  those  which  count  and  size  Individual  particles,  one  at  a time. 

The  former  type  Includes  the  well-established  transmlssometer  [7,8], 
which  measures  the  extinction  of  a light  beam  and  yields  a value  for  the 
Integrated  projected  area  of  the  particles  In  the  beam.  It  does  not  allow 
point  measurements,  and  requires  another  measure,  such  as  the  mass  loading, 
to  yield  a mean  diameter.  Other  techniques  have  been  described  [9-11]  which 
measure  the  angular  distribution  of  light  scattered  from  a large  number  of 
particles  simultaneously  present  In  a small  measurement  volume.  They  yield 
either  some  moments  of  the  distribution  [9,10]  or  require  the  angular  distri- 
bution to  be  unfolded  to  give  the  size  distribution  [11]. 
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Samp ling- type  optical  counter-sizers  for  gas  flows  have  been  commercially 
available  for  some  years.  The  sample  Is  diluted  as  necessary,  and  constrained 
to  flow  through  a small  jet  (~  1 mm  dla.)  which  Is  uniformly  Illuminated  with 
a collimated  white  light  beam.  The  scattered  light  pulses  are  collected, 
usually  In  forward  scatter,  detected  and  pulse-height  analyzed  to  give  a number- 
size  distribution.  Apart  from  a limitation  to  particles  5 10  ym  and  the  need 
for  sampling  and  sample  dilution,  the  main  problem  with  such  counters  [12,13] 
is  the  achievement  of  a response  which  is  insensitive  to  particle  refractive 
index. 

A novel  Interferometric  approach  to  in-situ  particle  counting-sizing  has 
been  suggested  by  Farmer  []4,15]  and  by  Durst  [16],  making  use  of  the  shape 
of  the  signals  scattered  by  single  particles  passing  through  the  fringes  in 
the  cross-over  volume  of  a dual  beam  laser  anemometer.  A simple  analysis  [14] 
suggests  that  the  signal  visibility  is  a monotonic  function  of  the  ratio 
particle  diameter/fringe  spacing,  in  a certain  range,  and  is  independent  of 
scattering  angle  and  refractive  index.  However,  more  detailed  analyses  [17] 
from  Mle  theory  and  experiments  [18]  Indicate  a complicated  dependence  on  these 
parameters,  which  raises  questions  concerning  the  method's  utility  except, 
perhaps.  In  a suitably  engineered  forward  scatter  geometry. 

Apart  from  this  interferometric  method,  there  appear  to  be  two  basic 
types  of  approach  to  devising  an  in  s^tu  particle  sizing  counter.  In  the  first, 
the  absolute  scattered  signal  amplitude  in  a single  collection-detection 
channel  is  employed,  as  in  commercial,  sampling  type  optical  counters.  In 
the  second,  use  is  made  of  the  relative  amplitudes  of  the  scattered  signals 
in  two  or  more  detector  channels  set  at  different  angles  [19]. 

In  this  paper  we  describe  the  development  of  an  in-situ,  forward  scatter 
laser  particle  sizing  counter  following  the  former  approach.  It  has  a demon- 
strated capability  on  cold,  small  scale,  low  velocity  flows  of  sizing  particles 

5 -3 

in  the  range  2 - 25  ym  at  concentrations  up  to  at  least  10  cm  . In  its 

present  form  the  Instrument  will  accommodate  flows  up  to  40  cm  in  width.  By 

adjusting  the  Instrumental  parameters  to  suit  the  particle  characteristics,  it 

should  be  possible  to  extend  the  range  downwards  to  ~ 0.5  ym  and  upwards  to 

6 *3 

i 50  ym,  and  handle  concentrations  to  10  cm*  . Moreover,  estimates  indicate 
that  the  Instrument  should  be  applicable  to  hot,  radiant  flows,  at  least  on 
laboratory  scale  systems. 
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There  are  two  major  factors  involved  In  the  practical  realization  of 
an  instrument  of  this  type.  The  first  concerns  the  control  of  the  size  and 
shape  of  the  effective  measurement  volume,  and  the  uniformity  of  response 
for  different  particle  trajectories  through  it;  this  is  discussed  in  Section  2. 
The  second  relates  to  the  attainment  of  a response  which  is,  as  near  as  possible 
monotonlc  with  particle  size,  and  which  has  minimal  sensitivity  to  particle 
refractive  index;  this  is  treated  in  Section  3.  As  discussed  later,  these 
two  factors  are  coupled  in  such  a way  that  the  optimum  trade-off  is  a function 
of  the  size  range  and  concentration  to  be  measured. 

From  the  discussion  of  Section  2 it  transpires  that  the  response  for 
3ifferent  trajectories  is  quite  non-uniform.  Hence,  as  experimental  tests 
show,  a flow  of  monodisnerse  particles  yields  a signal  amplitude-count  distri- 
bution with  a sharp  cut-off  at  some  maximum  signal  value,  and  a spread  to 
lower  values.  It  has  been  found  possible,  however,  to  use  a numerical  inversion 
scheme,  based  on  experimental  calibrations  with  monodlsperse  particles,  to 
unfold  this  distribution  of  signal  amplitudes  to  obtain  the  true  size.  This 
technique,  which  is  a vital  feature  of  the  present  scheme  is  developed  in 
Section  4. 

Section  5 details  the  results  of  calibration  tests  on  monodispersions 
and  polydispersions.  The  latter  demonstrate  the  instrument's  ability  to 
correctlv  resolve  the  structure  in  the  distribution. 

The  paper  is  concluded  with  a discussion  of  the  extension  of  the  technique 
to  a wider  size  range  and  its  application  to  hot  flow  systems. 
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2.0  MEASUREMENT  VOLUME  CONSIDERATIONS 

The  characteristics  of  the  measurement  volume  are  a crucial  factor  in 
the  design  of  an  in-situ  optical  sizing  counter  because  the  particle  trajec- 
tories are  not  controlled  as  they  are  in  a sampling-type  optical  counter. 

The  effective  measurement  volume  is  determined  in  part  by  the  intensity 
distribution  in  the  Illuminating  beam,  and  in  part  by  the  geometry  of  the 
collection  optics,  including  stops  and  apertures.  Ideally,  the  measurement 
volume  should  be  of  roughly  spherical  or  cubical  shape  with  a uniform  illumin- 
ation intensity,  and  a uniform  light  collection  efficiency  for  points  within  the 
volume.  Otherwise,  the  peak  scattered  Intensity  will  depend  on  trajectory, 
so  that  signals  from  large  particles  traversing  the  edge  of  the  volume  may  be 
comparable  with  those  from  small  particles  traversing  the  center. 

The  linear  dimensions  of  the  measurement  volume  should  be  several  (2  or 
3)  times  larger  than  that  of  the  largest  particles  to  be  measured,  so  that 
the  fraction  of  such  particles  which  follow  edge  trajectories  yielding  sub- 
standard signals  is  small.  For  a measurement  volume  satisfying  this  criterion, 
there  will  clearly  be  a limit  of  useful  operation  set  by  the  occurrence  of  a 
significant  fraction  of  coincidences  when  two  particles  are  simultaneously 
present  in  the  volume,  giving  rise  to  erroneous  signals.  This  condition, 
which  sets  a limit  to  the  maximum  mea.^urable  concentration  of  the  more  numer- 


ous small  particles,  is  analyzed  in  Appendix  I.  It  should  be  emphasized  that 
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the  capability  of  handling  concentrations  up  to  the  range  of  10  - 10  cm 


is  necessary  for  many  flow  systems  of  interest.  Roughly,  it  can  be  said 

3 

that  for  an  effective  measurement  volume  of  V cm  , the  maximum  concentration 

m 

that  can  be  measured  without  significant  Interference  due  to  coincidence  is 

-1  -3 

N ~ V cm  . 
max  m 

Concerning  the  distribution  of  intensity  in  the  Illuminating  beam, 
assuming  a TEM^^^ laser,  the  profile  is  Gaussian,  and  by  suitable  optics  can  be 
focused  to  a waist  at  the  measurement  station.  Using  cylindrical  coordinates 
centered  at  the  waist,  the  intensity  distribution  is 


I(r ,z) 


2 “ 2(r/w) 

IT  W l_  — 


(2.1) 
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where  P Is  the  beam  power  and  the  beam  radius  (to  1/e  Intensity)  is 


2 2 1/2 

w(z)  = w [1  + (Az/tt  w ) ] 
o o 


(2.2) 


Here  w is  the  beam  radius  at  the  waist  (z  = 0)  and  is  given  in  terms  of 
o 

the  far  field  half  angle  of  convergence,  0^,  by 


(2.3) 


Thus  a Gaussian  intensity  profile  (which  is  a r'^asonably  good  approximation 
to  the  ideal  rectangular  profile)  with  a waist  diameter  in  the  range  of 
Interest  (~  100  ym)  can  be  readily  created  by  suitable  beam  expansion  and 
focusing,  and  Involves  conveniently  small  beam  convergence /divergence  angles. 

It  may  be  noted  that  by  using  cylindrical  rather  than  spherical  lenses, 
ribbon  beams  with  Gaussian  profiles  having  distinct  waist  widths  in  orthogonal 
directions  (x,y)  can  also  be  created,  and  offer  certain  advantages. 

The  intensity  distribution  along  the  beam  (z)  axis  is  Lorentzian  and 
falls  off  much  more  slowly  than  with  radius.  The  Rayleigh  distance 
where  the  intensity  on  axis  falls  to  half  the  central  waist  Intensity,  is 


= TT  W /X 
O 


(2.4) 


For  waist  diameters  of  Interest  (2  w ~ 100  ym)  and  visible  lasers,  z„  ~ 1 cm, 

O R 

so  that  the  Illumination  volume  is  typically  very  long  compared  with  its  width. 

Clearly,  then,  it  is  necessary  to  limit  the  effective  length  of  the  measurement 

volume  in  the  z direction  by  suitable  design  of  the  collection  optics. 

The  most  direct  method  of  achieving  this  is  to  use  a lens  on  a collection 

axis  at  some  angle  6^  to  the  forward  direction  of  the  laser  beam,  and  image 

the  laser  waist  onto  a pinhole  aperture,  as  shown  in  Figure  1.  Clearly,  from 

the  point  of  view  of  limiting  the  length  of  the  measurement  volume,  a choice 

of  0 ~ 90“  is  most  effective.  On  the  other  hand,  as  discussed  in  Section  3 

c 

below,  to  satisfy  the  requirements  of  a monotonic  signal/size  response,  with 
least  sensitivity  to  refractive  index,  a coaxial  forward  scatter  geometry 
(0^  = 0)  is  most  desirable;  moreover  the  satisfaction  of  these  requirements 
becomes  more  difficult  as  the  scattering  angle  0^  is  increased. 


In  an  attempt  to  retain  a coaxial  forward  scatter  geometry,  and  to 
delimit  the  length  of  the  measurement  volume,  the  use  of  a central  stop  on 
the  collection  lens  was  investigated.  However,  a practical  solution  did  not 
appear  possible  with  this  configuration.  Consequently  It  was  found  necessary 
to  effect  a compromise  between  the  conflicting  requirements  of  limiting  the 
length  of  the  control  volume  and  achieving  a satisfactory  response,  by  using 
an  off-axis,  near-forward  scatter  geometry  with  0^  in  the  range  10  - 20®, 
as  shown  in  Figure  1. 

In  the  particular  system  investigated  the  illumination  is  provided  by 
a 2 mW  He-Ne  TEM^mode  laser,  which  Is  focused  by  cylindrical  lenses  to  a 
ribbon  beam  with  waist  widths  (to  1/e^  intensity)  of  approximately  2 w^  ~ 

100  ym  X 2 Wy  ~ 300  ym,  with  a waist  length  greater  than  1 cm.  The  collec- 
tion axis  is  in  the  x-z  plane  at  an  angle  0^  to  the  z-axis.  The  collection 

lens,  consisting  of  a back-to-back  pair  of  good  quality  f/2*8,  25  cm  focal 
length  lenses,  designed  for  infinite  conjugate  ratio,  images  the  center  of 
the  beam  waist  to  the  center  of  a pinhole  aperture  of  diameter  2 w^  ~ 100  ym 
diameter,  with  unity  magnification.  Thus  the  viewed  volume  is  approximately 
a cylinder  of  diameter  100  ym  intersecting  the  ribbon  beam  at  an  angle  0^, 
while  the  measurement  volume  is  given  approximately  by  the  intersection  of  the 
viewed  volume  and  the  illuminating  beam.  From  geometrical  considerations  the 
measurement  volume  can  be  characterized  roughly  as  a cylinder  with  slant  ends, 
of  diameter  2 w ~ 100  ym,  length  2 w /sin  0 ~ 1000  ym  for  angles  of 

rV  ^ ^ 3C  C 

interest,  and  volume  10  cm  . The  particle  flow  is  directed  in  the 

y-direction. 

The  foregoing  arguments  defining  the  measurement  volume  are,  of  course, 
only  approximate.  A rigorous  treatment  should  consider  the  angular  distribu- 
tion of  intensity  scattered  from  a particle  of  given  size  at  a given  position 
in  the  illuminating  beam,  and  thence  a calculation  made  of  the  intensity  dis- 
tribution in  the  plane  of  the  pinhole  produced  by  the  light  intercepted  by  the 
collection  lens,  to  obtain  the  detected  signal  as  a function  of  particle  size 
and  position.  Fortunately,  such  an  involved  calculation  is  circumvented  by 
the  calibration  procedures  described  in  Section  5.. 

In  any  event,  it  is  clear  that  the  use  of  forward  scatter  at  small  angles 
results  in  a distribution  of  intensity  in  the  measurement  volume  which  is 
quite  non-uniform,  although  the  use  of  a ribbon  laser  beam  reduces  the  non- 


uniformity  in  the  y direction.  Consequently  the  scattered  signal  is  trajectory 
dependent,  and  a flow  of  monodisperse  particles  yields  a signal  peak  amplitude 
count  distribution  having  a sharp  cut  off  at  some  maximum  signal  value,  corres- 
ponding to  particles  traversing  the  center  of  the  measurement  volume,  with  a 
spread  to  smaller  amplitudes. 

To  overcome  this  problem  of  nonuniformity  of  response  within  the  measure- 
ment volume,  a numerical  Inversion  scheme,  combined  with  a calibration  proce- 
dure has  been  devised,  as  discussed  in  Section  A,  to  unfold  the  distribution 
of  signal  amplitudes  and  yield  an  indicated  size  distribution  which  eliminates 
the  dependence  on  trajectory. 
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3.0  OPTIMAL  COLLECTION  GEOMETRY;  SCATTERING  THEORY 

Highly  desirable,  if  not  essential,  requirements  for  an  effective  optical 
sizing  counter  is  a monotonlcally  increasing  dependence  of  signal  amplitude  on 
particle  size,  and  Insensitivity  to  the  complex  refractive  index  n ■ n^  + ln2 
of  the  particle  (where  the  imaginary  pare  represents  absorption). 

To  find  an  optimal  geometry  for  satisfying  these  requirements  entails 
extensive  parametric  calculations  of  the  response  for  different  geometries  as 
a function  of  particle  size  and  refractive  index  from  Mle  scattering  theory. 

The  latter  gives  the  differential  scattering  amplitude  for  a spherical  particle 
in  an  incident  monochromatic  plane-polarized  wave,  as  a function  of  scattering 
angle  6,  polarization,  complex  refractive  index  n,  of  the  particle  (rela- 
tive to  the  medium)  and  the  normalized  size  parameter  a = ird/X,  where  d is 
the  diameter  and  X the  wavelength.  This  scattering  amplitude  must  then  be 
integrated  over  the  aperture  of  the  collection  optics  to  give  the  response. 

Unfortunately  the  Mle  scattering  function  is  an  extremely  complicated  and 
sensitive  function  of  the  parameters  6,  a,  n,  which  is  difficult  to  chara- 
cterize in  simple  terms,  although  fast  computer  codes  are  now  available  for 
its  evaluation.  However,  as  discussed  by  Hodklnson  [20],  in  the  limit  of 
large  particles,  a ^ 3/(nj^  - 1),  i.e.  d ~ 2X  ~ 1 ym,  a very  good  insight  can 
be  obtained  by  considering  the  separate  effects  of  diffraction,  refraction  and 
reflection.  Fifty  per  cent  of  the  total  scattered  power  is  attributable  to 
diffraction  and  is  distributed  in  the  well  known  and  easily  characterized 
Fraunhoffer  diffraction  pattern  of  a circular  disk  or  aperture,  which  is  the 
same  as  that  of  a spherical  particle  and  is  independent  of  polarization  and  of 
the  refractive  index. 

For  such  diffraction,  the  differential  scattering  efficiency,  defined  as 
the  flux  scattered  per  unit  solid  angle  divided  by  the  flux  incident  on  the 
particle  projected  area,  is 


2 r2  J 

4it  L a 


(a  sin  6) 
sin  0 


(3.1) 


where  is  the  Bessel  function  of  the  first  kind  and  order  one.  This  func- 

tion represents  the  familiar  central  lobe  surrounded  by  bright  rings 
of  decreasing  intensity.  The  first  zero  or  dark  ring  occurs  at 
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a sin  0 = 3.83.  For  a ~ 6,  and  X = 0.633  pm,  the  first  zero  occurs  for 
~ 44‘’/d(pm). 

The  Integrated  scattering  efficiency  is  given  by 


= 21: 


Idiff  sin  0 de 


= 1 - (a  sin  0)  - J?  (a  sin  0).  (3.2) 

o 1 


Approximately  85%  of  the  diffracted  flux  falls  in  the  central  lobe. 

This  result  leads  to  the  following  attractive  concept  for  particle  sizing; 
if  the  flux  is  collected  coaxially  by  a lens,  the  signal  will  be  monotonic 
with  size  and  independent  of  refractive  index.  Moreover,  for  a lens  aperture 


of  signal  will  be  closely  proportional  to  d for  d(pm)  > 

44°/0°  . For  instance,  for  0 ~ 10°  corresponding  to  F/2.8  collection 

in3x  in3x 


aperture,  this  will  hold  for  d IS  4 pm.  Of  course  there  is  also  an  upper 

2 

limit  to  the  size  which  yields  a signal  « d , in  view  of  the  need  to  stop 
off  the  center  of  the  lens  to  block  the  unscattered  beam. 

As  mentioned  in  Section  2,  the  use  of  such  a coaxial  forward  scatter 
geometry  was  investigated,  but  it  was  abandoned  in  view  of  the  difficulty  of 
limiting  the  axial  length  of  the  measurement  volume  for  reasonable  collection 
apertures.  Consequently,  to  obtain  a measurement  volume  small  enough  to  handle 
concentrations  of  practical  interest,  one  is  forced  to  accept  an  off-axis 
scattering  geometry. 

This  requires  computations  from  the  full  Mie  theory,  but  the  foregoing 
analytic  results  from  diffraction  theory  with  coaxial  collection  lead  one  to 
expect  that,  provided  one  keeps  the  central  scattering  angle  0^  small, 
then  the  response  will  not  depart  far  from  monotonicity  and  the  dependence  on 
refractive  index  will  be  weak.  One  also  anticipates  that,  in  general,  the 
use  of  a wide  collection  aperture  and  the  use  of  white  light  will  also  favor 
the  attainment  of  the  desired  requirements,  since  both  will  tend  to  average 
over  the  detailed  resonant  structure  characteristic  of  Mie  scattering. 

Lleberman  and  Allen  [21]  have  reported  experimental  and  theoretical 
results  for  the  response  of  a coaxial,  near  forward  scatter  configuration, 
corresponding  to  a commercial  sampling-type  optical  counter  employing  white 
light.  Scattered  light  is  collected  coaxially  between  angles  of  15.5°  and 
28°.  The  relative  response  as  a function  of  particle  diameter  for  latex 
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(n  = 1.6  + i.O)  and  glass  spheres  is  reproduced  In  Figure  2.  The  curve  shows 

2 

a generally  Increasing  dependence  on  d,  approximately  proportional  to  d , 
over  the  range  d = 0.3  ym  to  50  ym,  with  the  exception  of  a dip  or  resonance 
in  the  neighborhood  of  1 ym.  Specifically,  the  response  is  double  valued  for 
the  Interval  d = 0.65  ym  to  1.2  ym,  and  such  behavior  is  characteristic  of 
Instruments  of  this  type.  In  practice,  signals  from  particles  in  this  size 
range  are  lumped  into  a single  size  bin  in  the  pulse  height  analyzer.  For 
many  purposes  this  is  perfectly  adequate,  such  a 2:1  size  resolution  Increment 
being  comparable  with  that  achieved  by  cascade  irapactors  over  the  whole  size 
range. 

A systematic  calculation  of  the  response  curve  for  a coaxial  forward 
scatter  system,  as  a function  of  the  collection  aperture  between  scattering 
angles  and  02  (>  0j^)  was  made  by  Oeseburg  [22]  for  spherical  particles 

of  n = 1.49  + i.O,  and  for  size  parameter  up  to  a = 30  (corresponding  to 

d ~ 6 ym  for  X = 0.63  ym) . It  was  found  that  a resonant  dip  around  a ~ 5 

(d  --  1 ym)  was  a persistent  feature  over  wide  ranges  of  the  aperture  angles 
0j^  and  02,  and  that,  for  many  apertures,  additional  resonances  occurred  at 
several  larger  a values.  A measure  of  the  magnitude  of  a resonant  dip  can 

be  defined  by  R = 2 (ct^  - ag)/(aj^  + Qg) , where  ag  and  represent  the 

smaller  and  larger  size  parameters,  respectively,  which  give  double-valued 
responses.  Thus  R represents  the  size  resolution  Increment  for  which  signal 
amplitudes  must  be  lumped  into  one  size  bln,  and  has  a value  ~ 0.6  for  the 
case  discussed  above. 

Oeseburg  found  that  only  by  going  to  values  of  0^^  > 20°  combined  with 
very  large  values  of  02  - 120°,  could  the  R value  for  all  dips  be  signifi- 
cantly reduced.  The  optimum  aperture  for  the  chosen  refractive  index  occurred 
for  0j^  = 40°,  02  = 120°,  and  yielded  a response  with  no  dips  having  R 

values  greater  than  0.15. 

It  should  be  emphasized  that  while  such  very  large  apertures  are  feasible 
in  sampling-type  optical  counters  through  Che  use  of  ellipsoidal  mirror 
collectors  surrounding  the  measurement  volume,  they  are  quite  infeasible  for 
In-sltu  Instruments.  For  a useful  Instrumental  throw  between  the  measurement 
point  and  the  collecting  optics,  large  apertures  entail  expensive,  large  dia- 
meter lenses.  For  example,  the  collecting  lenses  used  in  the  present  vork. 
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RELATIVE  SCATTERING 


of  focal  length  25  cm  and  diameter  9 cm  (F/2,8)  allow  a maximum  difference 
(02  - 0j^)  ~ 20®  when  used  for  off-axis  collection. 

With  these  practical  restrictions  In  mind,  and  guided  by  the  general 
trends  of  the  results  of  references  21,  22  and  other  similar  published  work, 
response  calculations  have  been  carried  out  for  the  collection  geometry  Illus- 
trated In  Figure  3.  The  collection  lens  of  aperture  half  angle  0^  = tan  ^ 
(1/2F)  Is  centered  on  the  collection  axis  at  an  angle  of  0^.  Allowance  Is 
made  for  the  presence  of  a planar  mask  which  blocks  light  scattered  out  to 
an  angle  0^  > (0^  - 0l) , and  the  response  Is  calculated  by  Integrating  the 
Intensity  over  the  remaining  portion  of  the  lens,  shown  unshaded  In  Figure  3. 

The  calculation  of  the  response  function  proceeds  as  follows.  For 
spherical  particles,  van  de  Hulst  [3]  gives  the  basic  relation 


I = -fj  F’  (a,0,n)  • 
k r 


(3.3) 


where  (power /unit  area)  is  the  intensity  of  the  incident  plane  wave,  I 

is  the  Intensity  of  the  scattered  spherical  wave  at  radius  r from  the 
particle,  k = 2ir/X  is  the  wavenumber,  and  F'  Is  the  dimensionless  scatter- 
ing matrix. 

The  intensities  are  four-vectors  in  terms  of  the  Stoke 's  parameters 


I E 


r 

u 


(3.4) 


and 
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(3.5) 


For  an  unpolarlzed  incident  beam,  Eqn.  (3.3)  reduces  to 

^2  2 \ 2 I \ 


(3.6) 
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Figure  3.  Scattered  Light  Collection  Geometry 
for  Present  Optical  System 
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The  scattered  power  collected  in  an  aperture  of  area  A = r n,  where  Is 
the  solid  angle  for  collection  is  given  by 


F (a,n,n) 


(3.7) 


(3.8) 


where  the  dimensionless  response  function  is  defined  by 


(3.9) 


and  is  the  effective  scattering  cross-section  for  the  particular  collection 

2 

geometry,  normalized  by  X . 

The  matrix  elements  M^(a,6,n),  M2(a,9,n)  are  computed  by  a code  developr 

ed  by  Dave  [23].  The  response  function  F is  then  calculated  by  integrating 

(Mj^  + ^2)/^  over  the  collection  solid  angle. 

We  have  performed  calculations  on  65  sets  of  parameters,  and  although  not 

exhaustive,  the  results  from  these  calculations  indicate  general  trends.  The 

following  figures  (Figures  4-7)  are  representative  and  are  chosen  to  Illustrate 

the  typical  dependence  of  the  response  on  the  important  parameters.  To  reduce 

the  number  of  parameters  all  the  curves  shown  are  for  the  case  6 =9  i.e. 

me 

with  the  inner  half  of  the  lens  masked,  although  this  may  not  make  optimal  use 
of  the  lens  aperture.  Also,  an  unpolarlzed  beam  was  assumed,  though  the  choice 
of  a single  polarization  might  well  lead  to  an  Improved  response. 

Figures  4(a)-(c)  show  the  response  for  n = 1.6  + O.Oi  and  fixed  minimum 
scattering  angle  increasing  values  of  the  maximum  scatter- 
ing angle  (9  +9,)  i.e.  increasing  aperture.  Figures  5(a)-(d)  show  the 

c L 

response  for  the  same  refractive  index  and  fixed  aperture  (9  = 20“),  for 

1j 

increasing  values  of  scattering  angle  9^. 

All  of  these  curves  exhibit  a more  or  less  pronounced  resonance  (dip)  at 
a ~ 5,  together  with  a series  of  resonances  of  decreasing  magnitude  at  larger 
a.  This  form  is  insensitive  to  the  real  part  of  the  refractive  index,  n^^, 
especially  at  the  small  scattering  angles,  but  the  exact  positions  and  shapes 
of  the  resonances  do  depend  on  n^ . In  general  the  a values  of  the  resonances 
scale  as  (n^^  - 1)  . The  effect  of  increasing  the  aperture  to  Include  larger 

scattering  angles  (Figure  4)  is  to  reduce  the  magnitude  of  the  higher  resonances 
but  its  effect  on  the  first  resonance  at  a ~ 5 is  slight. 


15 
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Figure  5.  Variation  of  Response  Function  with 
for  6,  = 20°,  n,  = 1.6,  n„  = 0.0 
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In  view  of  the  persistence  of  the  resonances,  particularly  the  first, 
over  a wide  range  of  scattering  angles  and  collection  apertures  appropriate 
to  an  In-sltu,  forward  scatter  instrument,  one  should  examine  their  signifi- 
cance in  placing  limits  on  the  performance  of  such  an  instrument.  This  can  be 
discussed  in  terms  of  the  R values  of  the  resonances,  which  are  listed  for 
the  typical  case  of  Figure  4(b)  in  Table  3.1  below. 


Table  3.1.  Resonance  Widths  for  Figure  4b. 


“s 

a = (Og  + a^)/2 

d(ym)(for  X = 0.633  ym) 

R 

3 

7 

5 

1.0 

.80 

8 

13.5 

10.8 

2.2 

.51 

15 

19 

17 

3.5 

.24 

They  should  be  compared  with  the  R value  of  0.60  for  a commercial  sampling 
counter  as  discussed  by  Lleberman  [21],  and  with  the  value  of  0.67  (for  all 
particle  sizes)  typical  of  cascade  Impactors  having  a 2:1  incremental  size 
resolution. 

It  is  clear  that  an  instrument  with  an  acceptable  performance  with  respect 
to  resolution  can  be  realized  despite  the  presence  of  the  resonances.  The 
greatest  limitation  is  set  by  the  first  resonance,  which  is  only  relevant  if 
one  attempts  to  design  an  Instrument  to  work  through  the  value  ~ 5 
(d^  ~ 1.2  Um  for  X = 0.633  Um)  down  to  smaller  sizes.  The  presence  of  the 
higher  resonances  with  smaller  R values  should  present  no  significant  limit- 
ation to  measuring  the  larger  sizes  (d  > 2 m for  X = 0.633  ym) . For  many 
purposes,  where  high  resolution  in  the  range  1 ym  is  not  required,  a useful 
performance  may  well  be  obtained  even  through  the  region  of  the  first  resonance. 
Alternatively,  the  size  range  in  question  might  be  reduced  by  a factor  - 2 by 
the  use  of  a shorter  wavelength  e.g.  from  a He-Cd  laser,  which  would  also  have 
the  advantage  of  minimizing  Interference  from  thermal  background  in  hot  flows. 


As  noted  above,  for  the  range  of  scattering  angles  and  collection  apertures 
of  Interest,  the  dependence  of  the  response  on  the  real  part  of  the  refractive 
Index  over  the  limited  range  of  practical  interest  (n  ~ 1.3  - 1.7  say)  is 


^ 


weak.  However,  for  some  important  applications  e.g.  in  coal  combustion, 
there  can  be  an  appreciable  Imaginary  component  (i.e.  absorption).  For 
instance  a value  of  n2  ~ 0.1  corresponds  to  a highly  absorbing  material  such 
as  black  coal,  while  values  of  n2  “ 0.002  - 0.006  corresponds  to  measured 
values  [24]  for  the  sllicaceous  slag  in  a coal-fired  MHD  channel. 

Figures  6(a)-(d)  show  the  response  for  a fixed  collection  geometry 

0 =0  =4°,  0,  = 20“  and  n,  = 1.6,  for  increasing  values  of  n.,.  The 
mc  L 1 2 2_2 

difference  between  n2  = 0 and  n2  = 10  is  negligible,  but  for  n2  = 10  , 

the  slope  of  the  smooth,  large  a (a  > 20)  part  of  the  response  is  significantly 

reduced,  leaving  the  lower  a part,  including  the  resonances, virtually 

unchanged.  Further  increase  of  n2  to  10  ^ produces  a smoothing  of  the 

resonances  to  yield  a monotonic  response.  Figures  7(a)-(d)  parallel  those  of 

f Figures  6 except  that  a smaller  minimum  scattering  angle  ~ 

used.  The  main  difference  is  that  the  sensitivity  to  n2  of  the  slope  at  large 

[ a is  reduced.  Thus  for  large  particles,  a small  scattering  angle  should  be 

used  to  minimize  sensitivity  to  absorption.  i 

^ In  concluding  this  section,  certain  general  statements  can  be  made  as 

I follows.  For  practical  ranges  of  scattering  angle  and  aperture  appropriate  i 

j to  an  in  situ  Instrument,  an  adequately  monotonlc  response  function  can  be 

achieved  to  allow  useful  resolution  over  the  size  range  a > 5,  and  possibly 

to  smaller  sizes,  and  with  a relatively  weak  sensitivity  to  complex  refractive 

index  n.  However,  for  accurate  results  it  would  be  desirable  to  calibrate 

the  instrument  using  monodispersions  of  particles  with  the  refractive  index 

f of  interest.  In  general,  for  the  larger  sizes,  small  scattering  angles  should 

I be  used  to  minimize  sensitivity  to  n,  while  for  the  smaller  sizes,  larger 

I scattering  angles  should  be  used  to  minimize  the  measurement  volume  and  hence  j 

I raise  the  maximum  concentration  that  can  be  handled. 

i I 

i 1 
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Figure  6.  Variation  of  Response  Function  with  Complex  Refractive 
Index  n^  (Absorption)  for  = 1^,  0,  = 20*^,  n-  - 1.6 
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4.0  TRAJECTORY  DEPENDENCE  - INVERSION  TECHNIQUE 

As  discussed  in  Section  2.0,  In  devising  an  in-situ  partlcle-slzlng 
counter,  one  has  to  face  the  problem  that  the  peak  value  of  the  detected 
signal,  which  Is  used  to  size  the  particles.  Is  a function  not  only  of  particle 
size  but  also  of  Its  trajectory  through  the  measurement  volume.  The  principal 
factor  controlling  this  trajectory  dependence  Is  the  distribution  of  intensity 
of  Illumination  In  the  measurement  volume.  To  take  account  of  this  we  define 
I^(x,z)  as  the  peak  intensity  experienced  by  a particle  following  a y-directed 
motion,  as  a function  of  the  transverse  coordinates  (x,z)  of  its  trajectory. 
However,  the  detected  signal  peak  amplitude  is  also  determined  in  part  by  the 
collection  optics,  through  the  imaging  of  the  scattered  light  distribution 
incident  on  the  collection  lens  onto  the  pinhole  aperture.  This  Is  allowed 
for  by  defining  a dimensionless  transfer  function  H(x,z,a)  which,  in  general. 
Is  a function  of  the  particle  diameter  d (or  size  parameter).  In  practice, 
as  measurements  confirm,  the  function  H can  be  sensibly  independent  of  d, 
but  for  generality  it  is  included  by  allowing  for  the  combined  effects  of  non- 
uniform  Intensity  and  collection/imaging  by  defining  the  measurement  volume 
function 


J(x,z,a)  = I^(x,y)  H(x,z,a)  (4.1) 

The  detected  signal  peak  amplitude,  delivered  to  the  pulse  height  analyzer 
can  then  be  written 

.2 

A = G • — T(x,z,a)  F(a)  (4.2) 

4Tr'’ 

where  G is  a constant  representing  the  detector  sensitivity  and  electronic 
gain,  and  F(a)  la  the  response  function  defined  in  Section  3.0  for  particles 
of  size  parameter  a Illuminated  by  a uniform  plane  wave.  It  may  be  noted 
that,  especially  for  particles  large  compared  with  the  linear  dimensions  of 
the  measurement  volume,  the  particles  are  not  in  fact  illuminated  by  a uniform 
plane  wave.  However,  the  effect  of  this  is  taken  care  of  by  the  inclusion 
of  tne  size  dependence  in  J. 

Figure  8 represents  schematically  a quadrant  of  the  measurement  volume 
function  J(x,z,d).  It  is  shown  as  a monotonlcally  decreasing  function  of  x and 
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z measured  from  the  "center"  of  the  measurement  volume  where  It  has  its  maximum 
value  J^.  However,  the  following  analysis  is  valid  for  more  general  forms 
of  J. 

It  is  assumed  that  particles  have  an  equal  probability  of  passing  through 
any  element  of  the  cross-section,  and  that  all  particles  have  the  same  mean 
velocity  U,  Irrespective  of  size  and  trajectory. 

If  the  measurement  volume  function  J were  uniform,  with  value  J , 

m 

all  particles  of  the  same  size  parameter  (a,  say)  would  yield  the  same 

2 2 ^ 

signal  amplitude  = G(X  /477  ) J^FCOj^)  i.e.  the  quantity  A would  be  a unique 
function  of  a,  independent  of  x and  z.  In  practice,  if  a monodisperse  aerosol 
flows  through  the  measurement  volume,  and  a statistically  large  number  of 
signals  is  accumulated  in  the  pulse  height  analyzer,  a count  distribution  is 
obtained  having  a sharp  cut  off  at  the  maximum  signal  amplitude  A^,  and 
spreading  down  to  lower  amplitudes. 

Further,  if  a polydisperse  aerosol  is  used,  the  resulting  signal  ampli- 
tude distribution  consists  of  a superposition  of  such  monodisperse  distribu- 
tions. The  problem,  then,  is  to  devise  a technique  to  unfold  the  measured 
signal  amplitude  count  distribution  in  the  presence  of  a nonuniform  measurement 
volume  function  J,  to  yield  that  which  would  be  obtained  with  a uniform  J 
function,  and  hence  to  obtain  the  true  distribution  of  F and  from  this  the 
aerosol  size  distribution. 

Such  unfolding  or  inversion  techniques  have  been  applied  in  a number  of 
technical  areas,  notably  in  tomographic  X-ray  scanners  [25],  and  have  become 
technically  practical  since  the  advent  of  on-line  mini-computers.  Although 
a straightforward  mathematical  inversion  is  not  feasible,  given  the  experi- 
mental nature  of  the  problem,  efficient  solution  procedures  have  been  developed 
for  determining  a solution  which  has  an  error  magnitude  of  the  same  order  as 
that  of  the  measured  input  variables.  Further,  it  is  intuitively  plausible 
that  the  more  nearly  the  J function  approaches  the  ideal  rectangular  form, 
the  more  accurate  should  be  the  inversion  process.  Thus  there  is  an 
advantage  in  achieving  a form  for  J that  approaches  the  ideal  rectangular 
one  as  closely  as  possible. 

4 . 1 Formulation  of  Method 

The  problem  is  formulated  as  follows.  First,  it  is  convenient  to  normal- 
ize the  functions  J,  F and  A on  their  maximum  values.  Thus  we  define 
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J - (0  ^ J i 1)  . 


F - ^ (0  < F < 1)  , 

m 


(4.3) 

(4.4) 


where  F Is  the  maximum  response  function  value  corresponding  to  the 

HI 

largest  particle  In  the  aerosol  of  Interest,  or  the  largest  size  that  It  Is 
of  Interest  to  measure. 

The  maximum  signal  value,  corresponding  to  this  maximum  F value 
Is  then 


and  we  define 


A = G J F 

m ,2mm 

47r 


•f-  (0  < A < 1) 

A 

m 


(4.5) 


(4.6) 


Then  (4.2)  can  be  written  In  normalized  form 


J F 


(4.7) 


Now,  as  Is  appropriate  to  the  problem,  we  coarse-graln  the  functions  A,  F, 
treating  them  as  discrete,  rather  than  continuous  variables, and  let 


C(A^)  - ~ Signal  count  rate  for  normalized  signal  peak  amplitudes 


In  the  range  A^  to  (A^  + AA^) 


»(Fj) 


F N, 


concentration  of  particles  In  the  size  parameter  range 
yielding  normalized  response  function  values  In  the  range 
Fj  to  (Fj  + AFj) 


AS(Aj,Fj) 


= AS 


IJ 


cross-sectional  area  of  the  measurement  volume, 
normal  to  the  flow  direction,  which  yields 
normalized  signal  peak  amplitudes  In  the  range 
Aj^  to  (Ajj^  + AA^)  for  particles  having  normal- 
ized response  functions  In  the  range  F^  to 
(fj  > i?j). 


i 
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Then  the  count  rate  distribution  is  given  by  the  equation  set 


ni 


m 


E « 4 S N 
J“1  J ^ 


(4.8) 


m 

T “ 4 S , M, 
1-1  ”1  J 


which  can  be  written  more  compactly  In  matrix  form 


U A S • N 


(4.9) 


The  subscripts  1,  j are  written  in  order  of  Increasing  A and  F, 
starting  with  the  smallest  (utilized)  amplitude  channel  (1*1)  of  the  pulse 
height  analyzer,  and  the  smallest  response  function  class  (j“l)  capable  of 
recording  a signal  in  that  channel,  and  ending  with  the  largest  signal 
amplitude-class  (i=m)  and  largest  response  function  class  (j*m)  having  non- 
zero values.  The  discrete  classes  are  assumed  contiguous,  so  that 


^1+1  " ^^1  ^1^ 


and 


j+1 


(Fj  + AFj) 


(4.10) 


Hence  m is  the  number  of  (utilized)  channels  in  the  pulse  height  analyzer 
and  is  also  the  number  of  response  function  classes  into  which  the  particle 
distribution  is  resolved.  Here,  and  in  what  follows,  it  is  assumed  that 


F(a)  is  monotonlc,  at  least  when  coarse-grained  over  the  Intervals  AF 


j' 


The  signal  amplitude  classes  and  the  response  function  classes 


Fj  are  associated  on  a one  to  one  basis  through  the  choice 


Ai,  = J F, 
k m k 


or,  since  J 


(4.11) 
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since  a particle  yielding  a response  function  in  a given  class  Fj 
cannot  give  rise  to  a signal  amplitude  greater  than  A. , the  cross-section 


matrix  elements 
explicitly 


AS 


ij 


y 

are  all  zero  for  1 > j.  Thus  Eqn.  (4.9)  can  be  written 
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The  inversion  of  Eqn.  (4.12)  Is  written  symbolically 


N 


= i A S-^ 


U 


(4.13) 


Given  the  count  rate  distribution  of  signal  amplitudes  pulse 

height  analyzer  channels,  then  i^  the  AS  matrix  is  known,  Eqn.  (4.13)  can  be 
solved  to  yield  the  distribution  normalized  response  functions. 

If  the  response  function  F(a)  is  known,  either  from  computations  (when 
the  refractive  index  is  known),  or  from  a calibration  experiment  (see  Section 
5.0),  the  distribution  of  particle  size  parameters  Nj(aj)  and  hence  the  dis- 


tribution of  diameters  Nj(dj)  for  the  polydispersion  can  be  obtained.  The 
central  problem  is  the  determination  of  the  cross-section  matrix,  to  which 
we  shall  return  shortly. 

Up  to  this  point  the  functional  dependence  of  (and  hence  of  Fj^(k)) 

has  been  arbitrary.  Usually,  pulse-height  analyzers  offer  a choice  of  either 
a linear  amplitude  mode,  for  which  ^ ” constant  for  all  k, 

or  a logrlthmlc  mode  for  which 
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constant 


(4.14) 


i.e. 

= (1  + g)  (for  all  k) 


The  constant  g is  chosen  by  setting  the  gain  of  the  logarithmic  amplifier 
to  suit  the  relative  size  range  of  interest  and  the  resolution  i.e.  the  number 
of  channels  utilized. 

The  choice  of  the  logarithmic  mode  is  preferred  for  two  reasons.  First 

it  is  the  most  appropriate  mode  when  dealing  with  particle  distributions  with 

a large  relative  spread.  Second,  as  discussed  below,  in  the  case  that  J is 

independent  of  d,  a condition  that  can  effectively  be  realized  in  practice, 

a further  simplification  of  the  ^ matrix  results  which  reduces  the  number 

2 

of  independent  elements  from  m /2  to  m. 

With  a logarithmic  mode,  we  also  have  in  view  of  Eqn.  (4.11) 


— — M 

Fk 

i.e. 

K+1  " (1  + \ all  k) 


(4.15) 


4.2  Determination  of  the  Mi  Matrix 

Returning  to  the  question  of  determining  the  ^ matrix,  in  principle 
this  can  be  deduced  from  measurements  of  the  signal  amplitude  by  traversing 
a single  particle  of  known  size  dj  (e.g.  attached  to  a microscope  slide) 
through  the  measurement  volume  for  a range  of  x,  z and  d. 

Fortunately,  a simpler,  less  tedious,  and  more  direct  procedure  is 
available  through  the  use  of  an  aerosol  generator  capable  of  producing  mono- 
dispersions  of  controlled  size  and  concentration.  The  method  proceeds  as 
follows . 
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A monodisperse  particle  distribution  N(dj^)  with  a known  diameter  dj^ 
and  concentration  N Is  passed  through  the  measurement  volume  at  a known 
mean  velocity  U.  For  a specified  time  the  signals  are  accumulated  In  the 
pulse-height  analyzer  to  yield  the  count  rate  distribution  C(A^)  = C^.  Since 
N(dj)  = 0 for  all  j except  J = k,  the  system  of  equations  (23)  reduces 
to  a set  of  k equations  of  the  form 

C = U AS,,  N,  , (4.16) 

1 ik  k 

which  can  be  solved  directly  for  the  This  yields  the  k elements 

in  the  kth  column  of  the  matrix.  By  repeating  the  procedure  for  different 
size  classes  k,  successive  columns  of  the  ^ matrix  can  be  determined 
until,  by  m such  measurements  the  whole  matrix  is  known.  Incidentally, 

in  the  process,  the  cut-offs  of  the  count  distributions  at  Aj^  yield  an 
experimental  calibration  of  the  pulse  height  analyzer  channel  numbers  in  terms 
of  the  particle  diameter  (i.e.  F(d))  for  the  particular  refractive  index  employed. 

This  method  has,  in  fact,  been  used  and  is  satisfactory  though  time- 
consuming  if  m is  large.  In  the  process,  however,  it  was  realized  that  a 
much  simpler  technique  is  applicable  under  certain  conditions,  as  follows. 

When  the  amplitude  response  of  the  pulse-height  analyzer  is  logarithmic, 
then  if  the  transfer  function  for  collection  H,  and  hence  the  measurement 
function  J,  is  independent  of  d,  it  can  be  shown  (Appendix  II)  that  the 
AS  matrix  simplifies  so  that  all  the  elements  on  any  diagonal  are  equal.  It 
can  then  be  written 
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In  Eqn.  (4.17),  ASj^  Is  now  the  area  In  the  x-z  plane  where  the  J function, 

which  has  been  coarse-grained  into  logarithmic  intervals  with  AJ^/J^  = B,  has  the 

value  J„  (Figure  9).  Thus  AS  is  the  central  area  where  J has  its  maximum  value  i 

Jjji  = 1.  Also  is  the  peripheral  area  where  the  J function  has  its 

minimum  coarse-grained  value  = 1/(1  + 3)  , such  that  a particle  of  the 

largest  response  function  class  F generates  a signal  falling  in  the  smallest 

in 

signal  amplitude  class  A^. 

The  significance  of  this  simplification  is  that  now  all  the  matrix 
elements  can  be  determined  from  a single  calibration  using  monodisperse 
particles  of  the  largest  size  class  of  Interest,  from  the  equation  set 

C,  = U A S.  N , £ - l,m  (4.18) 

X x.  m 

The  calibration  procedure  is  then  particularly  simple  and  quick. 

Experimental  evidence  for  the  fact  that  the  J-functlon  is  sensibly 
independent  of  d,  at  least  for  the  conditions  of  the  experiment,  is 
adduced  in  Section  5.0,  and  hence  justifies  the  use  of  this  simple  calibration 
procedure. 

One  further  point  concerning  the  calibration  technique  should  oe  dis- 
cussed. To  determine  the  absolute  values  of  the  AS^  from  Eqn.  (4.18)  it  is 
necessary  to  know  both  the  velocity  U and  the  absolute  concentration  N^. 

In  practice,  when  using  the  monodisperse  particle  generator,  it  is  not  an  easy 
matter  to  determine  the  absolute  concentration  at  the  measurement  station. 

Even  if  the  particle  generation  rate  is  known,  it  is  difficult  to  allow  for 
loss  of  particles  by  deposition  In  the  feed  tube. 

If  only  the  relative  size-concentration  distribution  of  an  unknown 
aerosol  is  required,  then  relative  values  of  the  ASjj^  suffice,  and  these 
can  be  determined  without  knowing  either  U or  absolutely.  When  an 

absolute  calibration  is  required  the  concentration  N can  be  determined  by 

m 

an  extension  of  the  basic  calibration  technique  as  described  below. 

In  performing  the  calibration  procedure  based  on  the  use  of  Eqn.  (4.18) 
the  total  count  rate  in  channels  i = £ to  m is 

m m 

C.  = U V AS^  N (4.19) 
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Figure  9. 
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Schematic  of  AS^  Corresponding  to 
of  Figure  8 
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which  we  write 


= U s 


£,m 


m 


(4.20) 


C»  Is  determined  from  the  total  number  of  counts  In  channels  Jl  to  m 

Xj  yfn 

In  a known  time,  and  If  the  velocity  Is  measured,  we  can  determine  from 

(4.20)  provided  we  can  determine  the  area 


A S, 


(4.21) 


The  latter  Is  the  cross-sectional  area  of  the  measurement  volume  which  has 

J values  bounded  by  J = 1 and 
•'  m 


1 ^ 1 
(1  + e)™"^  (1  + 3)’' 


(4.22) 


where  r = (m  - H)  Is  the  difference  In  the  channel  numbers  over  which  the 
count  rate  Is  summed. 

The  area  ^ can  readily  be  determined  experimentally  by  traversing 
a small  pinhole  (e.g.  5 fjm  diameter) , much  smaller  than  the  detection  aperture 
pinhole,  through  the  measurement  volume.  The  extent  In  x and  z of  the 
cross-section  for  which  the  signal  falls  from  the  m^^  channel  (at  the  center) 
to  the  (m-r)^^  channel  Is  mapped  out,  and  hence  the  area  Is  found. 

In  this  way  the  absolute  concentration  Is  found  from  Eqn.  (4.20)  and 

then  the  absolute  values  of  the  AS^  are  found  from  Eqn.  (4.18).  The  use  of 
this  method  with  monodlsperslons  of  different  size,  to  check  experimentally 
that  J Is  Indeed  not  a function  of  d.  Is  discussed  In  Section  5.0. 


4.3  Implementation  of  Inversion  Technique 

Once  the  AS^^  matrix  elements  are  determined  from  such  a calibration, 
they  are  entered  In  the  matrix  Inversion  algorithm  of  a mini-computer.  Them, 
for  an  aerosol  under  Investigation,  the  count  rate  distribution  from 

various  channels  of  the  pulse  height  analyzer  are  automatically  entered  Into 
the  computer.  It  Is  useful,  however,  to  observe  the  data  on  the  Instrument’s 

visual  display  as  a check  on  correct  functioning.  ) 

1 
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4. 


Because  of  the  experimental  uncertainty  In  determining  and  AS, 
the  problem  N » (U  AS)  ^ C cannot  be  accurately  solved  by  straightforward 
matrix  inversion  techniques.  Instead,  the  solution  N is  derived  by 
searching  for  a non-negative  vector  N which  minimizes  the  residual  vector 
R = C - U ^ N.  A satisfactory  solution  to  this  problem  was  found  through 
use  of  a non-negative  least  squares  (NNLS)  solution  procedure  developed 
by  Lawson  and  Hanson  [26].  This  procedure  minimizes  error  propagation  and 
also  outputs  values  of  residuals  which  allows  one  to  judge  the  accuracy 

of  the  resulting  distribution. 
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5.0  EXPERIMENTAL  RESULTS 

The  objectives  of  this  section  are  two-fold: 

1.  Compare  experimental  results  for  the  scattering  response  function 
F(d)  with  the  calculations  of  Section  2. 

2.  Validate  the  accuracy  of  the  inversion  t;cheme  described  in 
Section  4. 

The  most  direct  method  of  validating  the  approach  described  up  to  this 
point  consists  of  measuring  well-characterized  monodispersions  of  known 
size  and  concentration  by  using  a Berglund-Liu  monodisperse  particle 
generator  [27].  Such  a device  is  claimed  to  generate  particles  within  ±1% 
of  the  stated  diameter.  Microscopic  examination  of  samples  collected  on 
glass  slides  confirm  monodlsperslty  to  within  the  ±5%  resolution  of  a 
microscope  graticule.  Spherical  droplets  are  generated  from  oleic  acid 
which  has  a refractive  index  of  n = 1.46  + O.Oi,  l.e.  negligible  absorption. 

For  a given  nozzle  diameter  (20  pm,  10  pm,  or  5 pm),  a 10:1  range  of 
particle  sizes  can  be  generated  by  varying  the  percentage  of  oleic  acid  in 
isopropyl  alcohol.  The  alcohol  evaporates  in  a long  residence  time  drying 

1/3 

chamber  leaving  oleic  acid  droplets  with  diameters  proportional  to  C , 
where  C is  the  oleic  acid  concentration.  The  present  experiments  used 
the  20  pm  nozzle  to  generate  particles  in  the  range  2.5  pm  to  25  pm. 

Tl'c  other  nozzles  (10  pm  and  5 pm)  can  be  used  to  generate  particles  down 
to  0.5  pm. 

Experimental  values  of  F(d)  are  derived  from  the  set  of  curves  shown 
in  Figure  ID.  Each  curve  corresponds  to  the  PHA  data  generated  for  a known 
particle  size.  The  relative  scattering  response  for  each  size  is  given  by 
the  threshold  channel  where  the  number  count  increases  sharply.  In 
practice,  there  exist  a few  doublets  and  triplets  of  the  primary  particle 
size  produced  by  the  particle  generator.  Thus  the  threshold  channel  is 
taken  as  that  channel  where  the  number  count  is  10%  of  the  peak  number  count. 
Comparison  of  these  experimental  results  and  the  calculations  of  F(d)  in 
Section  4 are  shown  in  Figure  U.  These  results  indicate  that  scattering 
theory  is  quite  adequate  for  the  range  of  conditions  investigated,  although 
further  work  is  necessary  to  confirm  the  extent  of  the  resonances  or  multi- 
valued regions  of  the  response  function. 
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Ignoring  the  resonances  and  taking  a mean  curve  fit  In  the  range  of 
2 - 30  pm  gives  an  approximate  scattering  function  dependence  on  particle 
size  as  F(d)  o d which  Is  similar  to  the  power  dependence  of  the  Royco 

PC  215  and  PC  245  particle  counters  where  the  exponent  ranges  from  1.5  to 
1,8  [12].  This  result  Is  to  be  expected  since  light  collection  geometries 
for  these  devices  are  similar  and  the  effect  of  a polychromatic  Illumination 
Is  only  one  of  smoothing  the  resonances.  Further  work  Is  continuing  to 
generate  particles  in  the  .5  - 2 pm  range  and  thus  provide  a check  on  the 
response  curve  for  small  particles. 

Proof  of  numerical  inversion  accuracy  requires  knowledge  of  calibration 

source  particle  number  density  along  with  particle  size.  The  number  of 

particles  produced  by  the  particle  generator  was  90,000/sec  and  the  the  flow 

conditions  used  in  our  experiments  gave  a maximum  number  density  of  108 
3 

partlcles/cm  . However,  a flexible  hose  was  used  to  convey  the  particles 
from  the  generator  to  the  measurement  volume  and  in  the  process  a 
significant  fraction  of  particles  was  lost  to  the  hose  wall.  In  addition 
it  Is  possible  that  the  local  number  density  traversing  the  measurement  volume 
may  not  correspond  to  the  bulk  mean  number  density  at  the  hose  exit,  even 
were  It  known.  Consequently  we  use  the  method  described  in  the  latter  part 
of  Appendix  II  to  determine  the  local  number  density.  It  Is  clear  that 
every  particle  above  the  minimum  size  threshold  traversing  the  measurement 
cross-section  will  be  counted,  as  long  as  coincidence  errors  do  not 

occur.  For  the  experimental  conditions  mentioned,  coincidence  errors  are 
negligible.  The  local  number  density  can  be  directly  evaluated  for  a 
monodlsperslon  by  the  formula  (Eqn.  II. 8). 


, ^k-r^k  _ 


where  ^ is  the  total  number  of  particle  counts  from  channel  (k-r)  to 
k where  r corresponds  to  from  Eqn.  II.  1, 

r * 2.3/8  for  8 < 0.2. 

In  general  S,  is  a function  of  particle  size  and  Table  5.1 

K“ir  y ic 

suinnarlzes  values  of  S.  . for  several  pinhole  diameters  and  two  values 
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Table  5.1.  Variation  of  Measurement  Volume  Cross-Section,  fc. 

With  P^nhole_Diameter , dj^,  and  Amplitude  Reference 
Level 


Pinhole  Diameter,  d, 
k 

^k-r, 

k 

j^^yHIIIIIIIII 

\-r'K  - 

\-r\  ■ -1 

5 pm 

.00525 

.00115 

12.5  pm 

— 

.0011 

25  pm 

.0045 

.00102 

The  values  of  ^ ^ same  for  the  5 pm  and  12.5  pm  pinholes 

while  r r <^®creases  by  11%  for  the  25  pm  pinhole.  Thus  for  particles 

< 25  pm  diameter  r r ~ ^r  sensibly  independent  of  particle  size. 

Given  that  is  reasonably  constant,  from  Appendix  II,  Equation  11.13, 

we  have  AS  ^./S,  = C^/C,  and  thus  the  normalized  cross  sections 

m-(k-i)  k-r,r  i k-r,r 

can  be  derived  directly  from  the  PHA  number  amplitude  spectrum  C^,  i = 1,  m. 

Figure  12  shows  a plot  of  normalized  cross  section  calibrations  for  five 

particle  diameters  ranging  from  6 pm  to  25  pm. 

A mean  curve  can  be  drawn  through  all  the  calibration  points  which  is 

within  ± 10%  of  all  values  for  the  five  particle  slzec,  indicating  that  the 

normalized  values  of  the  cross  section  matrix  are  sensibly  Independent  of 

particle  diameter  up  to  25  pm.  Coupling  this  result  with  AS,  s constant 

iC“i  y r 

up  to  25  pm  we  conclude  that  the  cross  section  elements  are  reasonably 
Independent  of  particle  size  up  to  25  pm  for  our  Instrument  geometry. 

We  know  that  the  beam  intensity  is  similar  for  other  overall  beam 

dimensions  (Eqn.  2.2)  and  thus  we  can  generalize  the  above  results  to  say 

that  the  particle  size  should  be  smaller  than  ' 35%  of  the  1/e  beam  width 

for  the  diagonal  elements  of  ^ to  be  identical. 

The  final  and  essential  test  of  the  inversion  scheme  and  its  accuracy 

is  to  obtain  data  on  an  aerosol  of  known  particle  size  distribution.  A 

mixed  aerosol  of  known  particle  sizes  was  obtained  by  sequentially  adding 

and  storing  the  data  from  four  monodlsperse  aerosols  generated  by  the 

Berglund-Llu  particle  generator.  This  accumulated  data  is  shown  in  Figure  13, 

which  represents  the  raw  PHA  data  generated  by  a mixture  of  6,  12,  17  and 
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Figure  12.  Non-Dimensional  Cross-Section  Calibration  Results 
for  Particle  Sizes  Ranging  from  2.9  - 25  um 


I2CX) 


Figure  13.  Pulse  Height  Analyzer  Data  for  a Mixture 
of  6,12,17,  and  25  um  Particles 


25  iim  particles.  Although  the  PHA  has  a total  of  100  amplitude  channels 


giving  6 >•  .06,  the  data  was  'compressed'  by  a factor  of  three  to  minimize 

computation  time.  The  'compressed'  channel  width  (P)  is  equal  to  0.21  which 

corresponds  to  a particle  diameter  resolution  of  Ad^^/dj^  = 0.13.  Application 

of  the  Inversion  scheme,  (Eqn.  4.13),  gives  the  results  shown  in  Figure  14. 
The  four  major  particle  size  peaks  occur  within  ±10%  of  the  known  particle 
diameters  although  a small  spurious  peak  occurs  at  channel  11.  The  peak  at 
channel  11  represents  a 10%  error  in  number  count  and  is  consistent  with  the 
accuracy  of  the  cross  section  calibration  curve  (Figure  12). 

An  approximate  analysis  of  error  propagation  in  the  inversion  scheme 
due  to  uncertainties  in  measurements  of  and  ^ is  presented  in 

Appendix  III.  The  results  of  the  error  analysis  show  that  the  upper  bound 
uncertainty  in  varies  from  ±13%  in  channel  9 to  ±26%  in  channel  20. 

Higher  counting  rates  for  practical  aerosol  flows  should  reduce  these 
uncertainties. 

-3 

It  should  be  remembered  that  the  generated  number  density  is  108  cm 
(an  upper  bound)  which  is  to  be  compared  to  the  computed  values  shown  in 
Figure  12.  The  actual  number  densities  at  the  nozzle  exit  of  the  particle 
generator  decrease  with  increasing  particle  size,  consistent  with  observed 
increases  in  particle  deposition  with  increasing  particle  size  [28]. 

The  computed  values  in  Figure  14  are  also  consistent  with  individual  particle 
size  calibration  tests,  Table  5.2,  where  the  number  density  was  calculated 
according  to  Eqn.  4.20. 

Table  5.2.  Comparison  of  Number  Density  Results  for 
Figure  12  and  Calibrations. 


Particle  Size  (pm) 

■BSraBI 

Calibration 

Fig.  12 
Results 

6 

73 

93 

12 

77 

66 

17 

63 

57 

25 

27 

31 

A 


41 


(g.ujo)  In  ‘AiiSNaa  daawnN 


42 


The  'calibration'  number  densities  also  generally  decrease  with  Increasing 
particle  size  although  exact  correspondence  of  the  number  density  with  the 
calibration  tests  was  not  found  because  the  deposition  rates  and  thus 
actual  number  density  at  the  nozzle  exit  are  very  sensitive  to  particle 
feed  line  curvature  which  was  not  carefully  controlled  In  different  tests. 


6.0  DISCUSSION  AND  CONCLUSIONS 

The  ^ situ  particle  sizing  counter  described  herein  satisfies  the 
principal  practical  requirements  of  such  an  Instrument  for  many  Important 
applications.  The  optical  access  requirements  are  modest  and  can  be 
readily  satisfied  for  laboratory  scale  flow  systems.  Particles  can  be  sized 
over  at  least  a 10:1  diameter  range  with  a resolution  Ad/d  - 0.15,  which  is 
adequate  for  most  aerosols  of  Interest.  The  Instrument  has  been  con- 
clusively demonstrated  as  capable  of  resolving  the  structure  of  a poly- 
dlsperslon  consisting  of  several  superimposed  monodlsperslons  in  the 
range  2.5-25  pm.  By  adjusting  the  Instrumental  geometry.  It  should  be 
straightforward  to  deal  with  sizes  ranging  up  to  -100  pm.  The  extension 
to  sizes  < 2 pm  appears  possible  but  requires  further  work  to  resolve 
the  limitations  set  by  the  first  resonance  in  the  scattering  function  which 
typically  occurs  In  the  range  of  1 pm.  Another  Important  feature  Is  that 
the  measurement  volume  can  be  made  small  enough  to  allow  measurements  with 
high  particle  concentrations,  in  the  range  of  interest  for  most  applications. 

The  key  feature  of  the  method  is  the  use  of  a computer-based 
Inversion  scheme  which  allows  one  to  unfold  the  effect  of  the  random 
trajectory  dependence  of  the  scattered  signals  together  with  a straight- 
forward calibration  technique  which  gives  the  necessary  Information  on  the 
scattering  function  F(d)  and  the  matrix  elemental  areas  required  for 

the  inversion. 

Although  the  experimental  work  was  based  on  the  use  of  non-absorbing 
spherical  droplets  supplied  from  a monodlsperse  particle  generator,  the 
method  should  also  be  applicable  to  absorbing  and  Irregular  particles. 

For  absorbing  spherical  particles  of  known  refractive  Index,  Mie 
scattering  computations  can  be  used  as  a guide,  but  If  the  refractive 
Index  Is  unknown,  or  the  particles  are  irregular,  the  calibration 
technique  can  always  be  used  provided  It  Is  possible  to  obtain 
monodesperslons  of  the  aerosol  of  Interest. 

Such  an  approach  has  been  described  recently  by  MarPle  and  Rubow[29,30] , 
with  application  to  Irregular  absorbing  particles  such  as  pulverized  coal, 
road  dust  and  other  aerosols.  They  used  a two-stage  Impactor  to  produce  a 
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sharp  upper  cutoff  in  particle  size  for  the  aerosol  of  Interest.  This  was 
then  fed  to  a sampling  type  optical  counter  (Royco  218)  to  give  a single 
calibration  point.  By  using  different  nozzles  in  the  Impactor  a series 
of  calibration  points  for  different  aerodynamic  diameters  was  obtained. 

Figure  15  shows  their  results  for  a variety  of  aerosols.  It  may  be  noted 
that  a mean  calibration  curve  could  be  used  for  all  of  these  aerosols 
without  introducing  a size  error  of  more  than  ±30%.  This  is  an  encouraging 
result  for  application  of  optical  sizing  to  various  Irregular  and  absorbing 
aerosols  of  practical  interest.  This  calibration  technique  using  an 
Impactor  is  directly  applicable  to  the  present  ^ situ  Instrument  for 
measurement  on  irregular  and  absorbing  particles. 

In  the  future  it  is  planned  to  explore  this  method  with  the  present 
in  situ  Instrument  to  determine  the  size  distributions  of  fly  ash  particles 
in  a cold  flow  electrostatic  precipitator  dust  tunnel. 

Another  potential  area  of  application  of  the  present  Instrument  is 
to  sizing  particles  in  hot  systems,  such  as  liquid  spray  combustors  and 
pulverized  coal  combustion  systems.  Two  potential  problems  may  be 
anticipated  in  hot  systems,  namely  Interference  from  background  thermal 
radiation  and  beam  and  image  wander  due  to  fluctuating  refractive  index 
Inhomogeneltles.  Estimates  of  the  effect  of  the  thermal  background  show  that 
this  should  be  negligible  for  sizes  down  to  0.5  pm  if  a laser  line  filter 
of  2o£  width  is  employed.  The  effect  of  beam  wander  is  more  difficult  to 
estimate  and  depends  on  the  amplitude  of  thermal  gradient  fluctuations 
expected  in  a particular  flow  system.  It  is  planned  to  study  such  effects 
by  making  measurements  in  a small  burner  system. 


ROYCO  218  OUTPUT  VOLTAGE  (volts) 
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APPENDIX  I 


SIGNAL  COINCIDENCE  EFFECTS 


For  a given  aerosol  size  distribution  and  instrumental  size  range,  there  is 
a limit  to  the  concentration  that  can  be  measured  set  by  the  occurrence  of 
interference.  The  latter  can  be  defined  as  occurring  when  the  signal  due  to 
one  particle  is  significantly  enhanced  by  signals  due  to  other  particles 
present  in  the  measurement  volume  "simultaneously",  i.e.,  within  the  resolving 
time  T of  the  electronics.  Clearly,  interference  is  most  likely  to  be  a 
problem  for  the  measurement  of  the  smallest  particles  in  the  distribution. 

Consider  a polydlsperslon  of  given  relative  size  distribution  (N^ (dj)/5^Nj) , 
and  suppose  the  instrument  is  set  up  so  that  the  maximum  size  class 
particle  will  register  signals  in  the  maximum  amplitude  channel  of  the 

pulse  height  analyzer,  while  the  lowest  amplitude  channel  A^^  just  registers 
signals  from  the  smallest  size  class  it  is  required  to  measure.  We 

wish  to  derive  criteria  for  how  large  the  total  concentration  can  be 

without  the  occurence  of  significant  Interference. 

For  simplicity,  we  take  this  to  be  the  condition  that  the  probability  of 
a signal  from  the  smallest  size  class  (1)  (strictly  the  smallest  scattering 
amplitude  class)  being  overlaid  by  one  or  more  additional  signals  of  equal  or 
larger  amplitude  is  less  than  some  critical  value  (say  ^ 20%).  Since  a 
simultaneous  signal  in  any  amplitude  class  i = 1 to  m will  cause  such  inter- 
ference, and  because  the  events  are  statistically  Independent,  the  required 
conditional  probability  is  just  the  probability  that  one  or  more  signals  is 
present  in  any  of  the  channels  during  the  resolving  time  x. 

Now  the  average  total  rate  of  signals  is  2^C. , and  the  probability  of  n 

i=l  i 

signals  being  present  simultaneously  is  given  by  the  Poisson  distribution 


P(n,x)  H (I.l) 

nl 

m 

where  x = y^C.  x 

i=l  ^ 

Hence  the  probability  of  one  or  more  signals  being  simultaneously  present 
is  1 - P(0,x)  = 1 - exp(-  x),and  the  required  probability  of  interference  is 
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? 


m 


Int 


1 - exp  -(  C t) . 
1=1 


(1.2) 


Since  we  require  that  is  small,  we  can  approximate  (1.2)  by 


^Int  “ ^1 


(1.3) 


Now,  from  Eqn.  4.8,  we  have 


y]c^  = U In,  AS  + N^  (AS  + AS  ,)  + - - 
i \ 1 m 2m  m-1 

+ N (AS  + + AS,)  ) 

mm  1 > 


m 


m 


-"E  Vj«) 

j=l  il=m-j+l  j=l 


»m 


(1.4) 


Here 


P.  ^ = Ut 
int 


j=l 


(m- j+1)  ,m 


(1.5) 


This  shows  that  the  larger  particles  in  the  distribution  have  a higher 
weight  for  causing  interference  because  there  are  more  channels  in  which  they 
are  capable  of  recording  interfering  signals.  Equivalently,  the  larger  particles 
experience  a larger  measurement  volume  j+1)  m^  ’ hence  have  a larger 

probability  of  detection.  It  is  also  clear  that  the  probability  of  interference 
depends  on  the  shape  of  the  distribution  as  well  as  the  total  concentration. 

In  Section  5 it  was  shown  that  for  a Gaussian  beam  profile  and  a logarithmic 
amplitude  mode  of  the  pulse  height  analyzer,  the  incremental  areas  AS^^^  are 
approximately  constant,  Independent  of  £,  so  that 


(m-j+l),m  = 


ASn  — j As 


(1.6) 


2’=m-j+l 
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and 


m 

''int  = ^ S J 


(1.7) 


This  can  be  expressed  in  terms  of  the  effective  measurement  volume  as  follows. 
Assuming  that  the  resolution  time  T is  equal  to  the  transit  time  U/L,  where 
the  length  L of  the  measurement-volume  in  the  flow  direction  is  set  by  the 
pinhole  diameter  in  the  collection  of  optics,  Eqn.  (1.7)  can  be  written 


P.  ^ ~ V (^)  N 

int  m \m/ 


j 


(1.8) 


where 


V = L S « Lm  AS 
m m 


(1.9) 

is  the  effective  (maximum)  measurement  volume.  Eqn. (I. 8)  can  be  written 

(1.9) 


P,  ^ ^ V 
int  m 


where 


m 


^ N. 

T pi  1 

is  the  total  concentration  over  the  size  range  being  measured,  and 

m 

"j 

1 = — 


(1.10) 


(I. 11) 


m 


j-i  ^ 

is  the  average  channel  number  determined  by  the  weight  of  the  distribution 

N.. 

1 

Equation  (I. 11)  shows  that  the  maximum  total  concentration  that  can  be 

measured  without  significant  Interference  is  determined  by  the  product  of  two 

factors,  and  (j/m)  which  we  discuss  in  turn. 

The  effective  measurement  volume  V depends  on  two  factors:  first,  the 

m 

optical  geometry  determined  by  the  waist  diameter  of  the  laser  beam  and  by  the 
F number,  collection  angle  0^  and  the  pinhole  diameter  of  the  collection 
optics;  second,  by  the  size  range  (‘^^lax^'^min^  measured,  as  it  determines 

the  range  of  scattering  amplitudes  hence  signal  amplitudes 
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A /A  to  be  accommodated.  If  S is  the  experimentally  determined  reference 
ml  _r 

area,  inside  the  contour  where  J = 1/e,  and  if  the  response  function  is 
approximated  by  F “ d^  (where  p ^ 1.65  for  our  experimental  conditions), 
then  we  may  write  for  the  effective  maximum  measurement  cross-section 


S 

ra 


P dn  (d  /d  , ) 
R max  min 


(1.12) 


In  practice,  the  optical  geometry  is  chosen  to  suit  the  maximum  particle 

diameter  to  be  measured.  For  Instance,  for  d ~ 25  pm,  the  geometry  used 

max 

to  obtain  the  results  reported  in  Section  5 resulted  in  a reference  area 

-4  2 2 

S ~ 4x10  cm  , l.e.  approximately  60  d . With  p = 1.65  and 
K mdx 

(d  /d  j ) = 10,  together  with  L = 100  pm,  we  find  an  effective  measurement 
max  min  ^ ^ 

volume  of  V ~ 1.5x10  cm  . 
m _ 

The  second  factor  (j/m)  is  determined  by  the  particle  distribution 

function.  Its  minimum  value  is  (j/m)  = (1/m)  when  the  distribution  is 

strongly  peaked  in  the  smallest- channel.  Thus  the  lower  bound  of  is 

V N_/m.  For  measurement  of  sizes  greater  than  the  mode  (maximum)  of  the 
m T 

concentration  distribution,  we  consider  distributions  of  the  form 


N(d)  = k d 

where  q is  positive.  This  yields 

.d^+Ad . 


N.  = 
3 


p - j 

1 N(d) 


k Ad. 
6d  ^ 

d.'l 

3 


For  F “ d^,  we  have 


hence 


Ad  « <6/p)  d . 

J ,^,-1 


In  the  special  case  q = 1,  we  have 


Nj  = k ($/p)  = constant, 


T = 1 V i = 

^ m 2-f  ^ 2 ’ 

j=l 


(1.13) 


(1.14) 


(1.15) 

(1.16) 


(1.17) 
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and  for  m » 1,  P,  N_/2. 

Int  mi  ^ 

In  practice,  most  distributions  decrease  more  rapidly  than  d . For 
instance  for  d greater  than  the  mode  of  the  distribution  by  mass,  q > 3. 
This  weights  j to  values  lower  than  m/2,  and  reduces  the  probability  of 
interference. 


In  general  for  N « d 


^(qFT) 

min 


^ I we  have 

k (g/p) 


(1.18) 


(d  /d  , ) 

max  min 


(j-i)(q-i)/(m-i) 


This  yields 


3 ~ 


(d  /d  , ) 
1 max  min 


(q-l)/(m-l) 


i(d  /d  . - (m+1) 

L max  min ' 


(d  /d  , ) 
max  min 


(q-1) 


J 


(1.19) 


For 


[(d  /d  . ) 
L max  min 


(d  /d  . ) 
max  min 


(q-l)/(m-l) 


« 1,  and 


(q-1)  > 


- m,  this  reduces  to 


T (gri.) 

^ ' (q-1)  Hn  (d  /d  , ) 
max  min 


(1.19) 


which  shows  that,  under  these  conditions,  the  factor  (j/m)  is  less  than 
(0.5)  and  decreases  as  the  distribution  falls  more  rapidly,  as  l/(q-l). 
Combining  (1.19)  with  (1.12)  then  leads  to 


P Si  — t 

int  (q-1) 


R T 


P Vr 


(1.20) 


Specifically,  for  (d  /d  , 

max  min 


(q-1) 

) = 10,  m = 20  and  for  q=l,2,3,4  we  have 


(j/m)  = 0.53,  0.33,  and  0.22  and  0.16  respectively.  Thus  the  factor 

(j/m)  decreases  rather  slowly  as  q Increases. 

In  particular,  for  the  Instrumental  conditions  reported  in  Section  5, 

-5  3 

for  which  V « 1.5  x 10  cm  , and  taking  Pj  _ = 0.2,  we  find  for  the 

™ 4 -3  4-3 

maximum  concentration  ~ 4 x 10  cm  for  q = 2 and  ^ 6 x 10  cm 

for  q = 3. 

For  d =25  tim,  and  a particle  specific  gravity  of  2,  these  dispersions 
tndx  ^ 

correspond  respectively  to  mass  concentrations  of  30  and  20  grams/m  , 
which  are  rather  higher  than  the  fly  ash  loadings  at  the  exhaust  of 
pulverized  coal  combustors. 
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APPENDIX  II 


CONDITIONS  FOR  IDENTITY  OF  DIAGONAL  ELEMENTS  OF  AS 


With  the  choice  of  a logarithmic  response  on  the  pulse-height  analyzer 
we  have  from  Eqs  (4.10),  (4.11)  and  (4.14)  the  relations 


Vi  - (1  + 6)  \ - Vi  ' « * ’v 


(II. 1) 


valid  for  all  k.  To  be  consistent  with  Eqn.  4.7,  J is  coarse  grained 
Into  levels  (1  - H - m)  such  that 


= (J  + AJ  ) = (1  + 6)  J 


(11.2) 


For  specified  3 and  given  there  Is  an  associated  cross-sectional  area 
AS^  which  In  general  Is  also  a function  of  the  particle  diameter  d.  However, 
for  the  case  where  AS^  Is  Independent  of  d,  a simpler  form  for  ^ results. 

In  the  equation  set  (4.12),  the  contribution  to  from  Nj  (j  > 1)  Is 
U AS^j  Nj . Such  contributions  can  be  assumed  to  be  generated  In  the  region 
of  the  measurement  volume  with  value  J having  area  AS, , so  that 


^1  ■ 

But,  In  view  of  relations  (II. 1)  and  (II. 2)  we  also  have 


(II. 3) 


A = J F 
*1+1  I *^j+l 


Hence,  It  follows  that 


AS 


ij 


AS 


1+1,  j+1 


AS, 


i.e.,  all  the  elements  along  any  given  diagonal  In  the  AS 


iJ 


(II. 5) 

matrix  are  equal. 


Furthermore,  since  It  Is  clear  that  I Is  constant  If  (J-1)  Is  constant,  and 
In  view  of  (II. 3),  then  the  relation  between  the  Indexing  schemes  Is 
H = m -(j-1) , I.e. , 


AS. 


AS 


(II. 6) 


'Ij  “ m-(j-l) 

Thus  we  have  shown  that  when  J Is  not  a function  of  d and  Intervals  of 
AAji^,  AFj  and  AJ^^  are  chosen  to  be  logarithmic  (Eqns  II.  1 and  II. 2),  the 
general  form  of  the  AS  matrix  appearing  In  Eq  (4.12)  reduces  to  that  of 
Eq  (4.17). 

Physically  this  means  that  for  all  pairs  of  (l,j)  values  such  that 
(i-j)  = (m-J,),  signals  in  the  amplitude  class  arising  from  particles 
yielding  response  functions  In  the  class  F^  are  generated  In  exactly  the  same 
part  of  the  measurement  volume  AS^^  corresponding  to  the  value  J » J^. 

The  question  arises  of  what  physical  conditions  have  to  be  satisfied  In 

order  that  J should  be  Independent  of  d,  and  how  can  one  check  experimentally 

that  the  condition  Is  fulfilled.  Clearly,  one  necessary  condition  Is  that  the 

largest  particle  to  be  measured  (In  the  class  d^  corresponding  to  the  response 

function  class  F ) should  be  reasonably  small  compared  to  the  dimensions  of 
m 

the  measurement  volume.  In  particular  the  radius  w^  of  the  laser  beam.  Other- 
wise, such  particles  will  not  experience  an  Incident  wave  that  Is  approximately 
a uniform  one,  as  Is  required  by  the  Mle  theory.  Most  likely  there  are  other 
conditions  related  to  the  Imaging  of  the  particle  onto  the  collection  aperture 
which  must  be  satisfied  In  order  that  J should  be  Independent  of  d.  However, 
these  are  difficult  to  specify,  since  they  probably  involve  the  quality  of 
the  collection  lens  with  respect  to  abberatlons. 

In  view  of  such  uncertainties  In  specifying  the  necessary  and  sufficient 
conditions  that  J should  be  Independent  of  d,  we  have  taken  a pragmatic 
approach,  and  devised  an  experimental  test  which  confirmed,  at  least  for  the 
experimental  conditions  employed,  that  J was  Indeed  sensibly  Independent  of  d. 
This  experimental  test  can  be  explained  as  follows.  At  the  end  of  Sec.  4.2, 
an  extension  of  the  basic  calibration  technique  was  described  that  allows  the 


i 
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absolute  concentration  N and  hence  the  absolute  values  of  the  AS.  matrix, 

m I 

to  be  determined,  using  the  largest  size  monodlsperse  particles  of  Interest. 
The  procedure  Involved  summing  the  count  rates  from  channels  1 = m down 

to  1 = (m-r) , l.e.,  over  the  top  (r+1)  channels,  and  using  the  relation  (see 
Eqn.  4.19) 


C = U S N 

m-r,r  m-r,m  m 


(II. 7) 


The  value  of  S , the  cross-sectional  area  of  the  measurement  volume  over 
m— r ,m 

which  J lies  between  J = 1 and  J = l/(l+6)''  is  found  by  traversing 

m m— r 

a small  pinhole  through  the  measurement  volume.  The  test  of  the  Independence 
of  J on  d consists  of  repeating  this  procedure  for  different,  small  sizes 
of  monodlsperse  particles  (Table  5.1). 

For  particles  of  size  class  d^,  and  concentration  Nj^,  the  sum  of  the 
count  rates  In  the  (r+1)  channels  from  1 = (k-r)  to  1 = k Is  given  by 


Vr.k  = U Vr.k  \ 


(II. 8) 


provided  that  J is  not  a function  of  d.  Under  the  same  condition,  If  r 
Is  held  constant  in  the  two  measurements. 


S 


k-r  ,r 


S 

m-r,r 


(II. 9) 


so  that  the  same  area,  determined  using  the  pinhole  traversing  technique,  can 

be  used  In  Eqn.  11.8  as  In  Eqn.  11.7.  Hence  the  absolute  value  of  Is 

determined  in  the  same  fashion  as  N . 

m 

Now  the  measurement  of  N allowed  the  absolute  values  of  the  AS, 

m i, 

determined,  for  = 1 to  m from  (Eqn.  4.16) 


(II. 10) 


Likewise,  the  measurement  of  allows  an  independent  set  of  absolute 

values  of  the  AS^^  to  be  determined  for  £ = (m-k)  to  m from  the 
equation  set 

= U ASj^  Nj^  (II.  11) 

where  £ = m-(k-l),  and  1 ranges  from  1 to  k.  It  follows  from 
Eqns.  II. 8 and  II. 10  that 


I 
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i = l,k  . 


(11.13) 


AS  C. 

m-(k-i)  _ 1 

S,  , C,  , 

k-r,k  k-r,k 


Thus  values  of  AS  „ normalized  on  a fixed  reference  area  S.  ^ _ 

can  be  determined  without  knowledge  of  the  number  concentration  Nj^.  This 
equation  is  valid  even  if  J has  a particle  size  dependence,  as  long  as  the 
reference  area  S,  is  determined  for  each  particle  size  d . For  the 

K“ir  j IT  iv 

case  where  J is  independent  of  d^,  then  k ~ ^r  ^k’ 


1 


APPENDIX  III  ERROR  ANALYSIS  OF  INVERSION  TECHNIQUE 


There  are  two  sources  of  errors  in  computing  values  of  particle  number 
density,  N^,  namely  random  counting  errors  In  and  errors  In  specifying 

correct  values  for  the  elements  of  Counting  errors,  which  derive  from 

a finite  sampling  time,  can  be  readily  evaluated  according  to  Poisson 
statistics.  We  have  assumed  that  all  particles  and  particle  size  classes 
have  an  equal  probability  or  passing  through  any  point  of  the  measurement 
cross-section.  This  assumption  Is  true  only  as  C^  approaches  Infinity. 

For  a 95%  confidence  level  Poisson  statistics  give 


where 


= random  error  In  C^ 


(III.l) 


C^  = measured  number  count 

According  to  the  above  relation  one  should  acquire  a minimum  of  1000 
counts  in  each  size  range  channel  to  achieve  a t6%  error  in  C^. 

The  effects  of  these  counting  errors  on  the  uncertainty  In  N^  cannot 
be  quantified  without  making  some  approximations  for  the  form  of  AS.  For  a 
Gaussian  illumination  beam  It  turns  out  that  the  Individual  cross-section 
elements  are  all  nearly  equal.  The  following  analysis  assumes  that  they 
are  identical  and  known.  For  = AS  for  all  1,  j the  equation  set  4.17 

becomes 

c,  = (Nt+n,+  . . . . + n ) u as 

C-  = N,+  . . . . + N ) U AS 

2 2 m 

• . • (III. 2) 


C = 
m 


(N  ) U AS 
m 


To  simplify  the  problem  we  choose  a velocity 
there  were  no  uncertainties  In  the  Individual 


U such  that  U AS  •»  1.  If 
A S^j  then  we  would  obtain 


”i  ° ^1  "‘^l+l’  ^ 


(III. 3) 
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The  uncertainty  In  N would  be  given  by 


(III. 4) 


or 


AN^ 

nT 


>JT 


1 


If 


< 1000  then 


’i  '^'^i  ^'i+l 

[aNj^I^N^  i 0.09  from  random  counting  errors. 

The  potential  errors  in  the  computed  number  distribution  N due  to  uncer- 
tainty in  the  elements  of  ^ require  numerical  calculations  using  perturbed  valu 
of  As  shown  in  Figure  12  an  average  curve  with  a deviation  of  ±10% 


characterizes  AS 


m-(k-l) ' 


Although  values  of  the  true  AS  elements  would  be 


expected  to  vary  systematically  with  Increasing  particle  size,  we  will  assume 
that  the  measurement  errors  are  random  within  the  range  of  the  five  sets  of 
experimental  results.  This  approach  will  provide  an  upper  bound  for  the  error 
estimate. 

For  numerical  analysis  of  the  effect  of  uncertainties  in  AS,  artiflcally 
perturbed  values  of  ^ were  generated  by  using  the  following  equation, 

~ (III. 5) 


where 
deviation  from 


AS^j  = AS^j  (0.9  + 2eRN) 

AS^j  is  a mean  experimental  value,  e is  the  maximum  fractional 

AS...  ^ind  RN  is  a random  number  between  zero  and  one.  Three  sets 


Ij’ 


of  ^ were  generated,  and  using  a given  count  distribution  C,  three  number 

distribution  solutions  N,  shown  in  Table  III.l  were  obtained.  The  overall 

standard  deviation  between  the  average  of  N for  the  three  perturbed  solutions 

— -3 

and  the  N for  4S^j  = ^S^^  is  ±5  cm  . Although  greater  deviations  from  the 

average  occur  for  the  perturbed  solutions  (e.g.  channel  13),  it  should  be 

noted  that  these  large  deviations  are  compensated  by  a corresponding  deviation 

of  the  opposite  sign  in  an  adjacent  channel.  In  other  words,  size  resolution 

and  number  resolution  are  coupled.  Forcing  better  size  resolution  produces 

% 

greater  errors  in  the  number  density  for  each  size  range  and  vice  versa.  The 
upper  bound  for  combined  uncertainty  from  random  counting  errors  (c.e.)  and 
uncertainty  in  ^ is  given  by 
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F 


Table  III.l.  Comparison  of  for  Perturbed  and 

Average  Values  of  AS. 


' ^ ' CE  ' ^ ' AS 


(III. 6) 


For  Che  results  of  Figure  14,  the  greatest  uncertainty  occurs  In  channel  20 
and  Is 


- ±26%. 


For  larger  number  density  flows  (as  occur  In  practice)  or  for  longer  sampling 
times,  counting  errors  can  be  significantly  reduced.  Reduced  counting  errors 
will  also  reduce  the  uncertainty  In  AS. 
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